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Abstract Cal span Corporation 

Given an unknown system with n inputs and a outputs, one possible 
node ling approach is to compute estimates of transfer functions using the 
mathematical methods of tfiener-Kopf (i.e., the t toady state Hainan as tins tors). 
Unfortunately, all the computational methods f r r implementing the Wioner-Hopf 
approach require (or assume) an analytical expression for the determinant of 
the spectral matrix. These assumptions, in turn, pre-ordain the final form of 
the Wiener estimates and hence the approach tends to lose the "Froe-Minimltatlofl 
methodology which makes it attractive. Clearly, what is needed is an algorithm 
which operates on the finite length input and output signals in a maimer which 
does not bias the results through the use of possibly unwarranted assumptions. 
7he procedure set forth in this paper furniskvs such an algorithm for tha 
special case where n finite-length realisations, for each of the n input signals 
are available. 

1 . Introduction 

The tfiener-Hopf method leads to tha formalisation of a spectral matrix 
which must be factored into the product of two matrices, one of which is analyti 
in the right half of the complex frequency plane; the other analytic in the left 
half plan*. Unfortunately, all methods for accomplishing this factorisation, at 
one particular step or another, require (or assume) an analytical expression for 
the determinant of the spectral matrix. This assumption, in turn, pre-ordain* 
the final form of the Wiener estimates. 

For example, if it is assumed that the determinant of the spectral 
matrix is a rational polynomial in 3 , then the estimators are forced to be 

rational polynomials in S also. Clearly what is needed is a computational 
algorithm which operates on the finite -length input and output signals and 
furnishes optimal estimates without prejudicing *he results through the use 
cf possibly unwarranted assumptions. The procedure set forth in this paper 
furnishes such an algorithm for the special case where the unknown system has 
n inputs and n finite length realization* for each of the n input signals are 
available. 

2. Theory 

In the interests of brevity, the results will be developed for the 
case where the unknown system has the two -input and one-output configuration 
shown in Figure 1. After this, the extension to the case of n inputs and m 
outputs will be obvious. 
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The Nlenor-Mopf minimi tat ion procedure loads, for the ideal case of 
an infinite number of infinite length realisations of tha signals, to Equation (1): 
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( 9 ) 


The details of solving this equation analytical ly ara discussed, for exeaple, 
in Rsfsrsncs 1 and 2. 

For the case whore only a finite nuaber of fixed length records are 
available, the various power spectra in Equation (1) are reoiaced by their 
estimates. For exaaple, the estimate of pi , cell it p , is 


K * £, < s > 


Suppose now that only two real iaat ions of each signal are available. 
Equation (1) becoaes, using j? to denote , etc.: 



Because we have chosen to work directly in the frequency doaaln with the 
transform of the various signals, the factorisation of the two input case becoaes 
trivial when only two realisations of the inputs are available. Factor Equation 
(6) as 
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Coapute the inverses of the 0 and 0 § given in equation (7) and substitute 
into equation (4) 




Thus one need only work with the transform# of the sign*]* when only 
two realizations are available. Note that no assumptions were necessary - the 
Wiener- Hopf theory guarantees that the H's computed according to Equation (9) 
are the best mean square fits possible under the stated conditions. 

The extent ion to the n input case is now obvious but will require e 
double subscript notation (refer to Figure 2). 


i (an m vector) 


FIGURE 2. DOUBLE SUBSCRIPT..' INPUTS 



The factorization is now 0 where 
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FIGURE 3. BLOCK DIAGRAM OF DIGITAL SIMULATION 

In Figure 3, ** and are wh*te noiae sources fed through two different 

filters to produce the input signals X and Y. The estimates of ft and '/**, were 
then computed using Equation (9). The results are shown in Figure 4 for three 
different values of the ratio of *//*j . At can be seen, the estimate of ft 

la excellent when the n, noise dominates while the estimate of '/**f is excellent 
when dominates. Somewhere in between, poorer estimates of both t/ A and 
are obtained. * 

The experiment also provided an opportunity to check out the effect 
of operating on the data with threa diffartnt time windows. These are denoted 
at the square window ( n }, the triangular data window ( a ) and the triangular 
window convolved with itself (AAA ). The resulting magnitude plots for 
f /s (Xf/Xg* ih ), shown in Figure 5, demonstrates the dramatic effect tV. 
time window has on the variability of the estimates. 


Conclusion 


Given an n input system and n expeiimental records for each input, it 
has been shown that the optimum Wiener filter is easily computed using only the 
experimental data. That la, no assumptions concerning the analytical structure 
of the spectral matrix are necessary. 
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